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ABSTRACT 


The numerical model of Williams et al. (1992) is used to study frontogenesis 
from unbalanced initial conditions. The dependent variables are assumed to be 
independent of y. The hydrostatic Boussinesq primitive equations are used with no 
diffusion of heat or momentum. The model is bounded at the top and bottom by rigid 
planes. Periodic boundary conditions are used in the horizontal. The lateral boundaries 
are placed far enough from the imbalance region to avoid wave reflection. The 
atmosphere is assumed to have constant vertical temperature stratification. 

The initial imbalance is obtained by allowing a horizontal temperature gradient 
to exist while the initial wind is zero. In a stably stratified atmosphere, gravity waves 
are excited and propagate away from the imbalance region, provided no reflection 
occurs in the lateral boundaries. Therefore, the atmosphere tends toward a geostrophic 
balance away from vertical boundaries. Near these boundaries, the temperature 
gradient oscillates or it collapses into a front, depending on the initial Rossby (Ro) 
and Froude (F) numbers. A relationship between Ro and F is established which 
separates situations where a front may or may not form. Numerical solutions show the 
formation of a front within a finite period of time that tilts toward the cold air. 
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I. INTRODUCTION 


Frontogenesis can be described as the process of formation of small-scale motion 
field from an initial larger scale state. Baroclinic instability is the main source of 
disturbance energy, that can eventually produce the energy cascade from the large to the 
small scale atmospheric states (Williams 1967; Hoskins and Bretherton 1972). Many 
different mechanisms can work to produce or modify a front. Numerical experiments 
(Williams et al. 1992) show the weakening of a front approaching a topographic ridge and 
reintensifying as it goes downslope. Turbulent horizontal and vertical diffusion can play 
an important role in opposing frontogenetical effects and in modifying the frontal 
structure (Williams 1974; Keyser and Anthes 1982). Background shearing or stretching 
deformation wind fields can generate large temperature gradients and are equally 
important in frontogenetical processes (Stone 1966; Hoskins and Bretherton 1972). 
Differential heating can be either frontogenetical or frontolytic, depending on whether it 
enhances or weakens the temperature gradient. 

Balanced frontogenesis has mass and momentum fields in near thermal wind 
balance and relative vorticity is larger than horizontal divergence. Williams (1972) used 
balanced initial conditions to show frontal formation with linear and nonlinear primitive 
equation models. The latter is more realistic in the sense that a discontinuity occurs 
within a finite period of time and the front tilts toward the cold air. This is due to the 
temperature advection by the divergent part of the wind, not present on the linear case. 
Hoskins and Bretherton (1972) suggest that relative vorticity will remain larger than 


1 



horizontal divergence during this process, in the absence of friction, even when the 
Rossby number is greater than one. 

The frontogenesis is unbalanced when the mass and momentum fields are not in 
semi-geostrophic balance. This situation is usually associated with large near-surface 
differential heating. A large temperature gradient can cause a non-geostrophic horizontal 
pressure gradient that can lead to frontal formation. The imbalance can originate in strong 
baroclinic zones associated with the sea-breeze circulation and in regions where partial 
cloud cover can cause large temperature gradients. Clear skies ahead of and extensive 
cloudiness behind a cold front can produce an imbalance during the morning heating. 
This imbalance can lead to a stronger front than the balanced one, which is forced by the 
large-scale flow. 

Blumen and Williams (1996, hereafter BW96) investigate unbalanced 
frontogenesis for a zero potential vorticity case. In one of their initial conditions, the 
temperature has an error function distribution and the fluid is at rest. A heat source can be 
considered to generate the imbalance but no somce exists after the initial time. The 
equations use the Boussinesq approximation and friction is neglected during the process 
of frontal formation. They find an inertial oscillation in the horizontal velocity 
components, and frontogenesis when the Rossby number is sufficiently large. 

The purpose of this study is to extend the work of BW96 for the case of a non¬ 
zero constant potential vorticity, with similar initial conditions. When the atmosphere is 
vertically stratified, gravity waves are allowed to propagate away from the initial 
imbalance and the frontogenesis process can be significantly modified. The main 



objective of this research is to determine a curve in Rossby number - Froude number 
space which separates frontogenesis from non-ffontogenesis cases. Various analyses and 
displays will be used in order to determine whether a frontal discontinuity has formed. 
The boundary conditions are the same as used by BW96, and heat sources and friction are 
neglected. Williams and Kurth (1976) also considered this problem, but with different 
initial conditions. 
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n. MODEL FORMULATION 


A. BASIC EQUATIONS 

The study of surface frontogenesis in this work involves just the lower portion of 
the troposphere where the density variation is small. To simplify the problem of modeling 
the frontal formation, the atmosphere will be considered incompressible, inviscid, and in 
hydrostatic balance. The equations developed here follow Williams et al. (1992) adapted 
to a flat topography. The vertical axis is normalized with the height of the upper 
boundary, H, and defined as 


Using equation (1) to transform the Boussinesq equations (Ogura and Phillips 
1962) we get 


dW ^dw 

—+V.VV+C^ = -V0-^xV+F, 
ot oQ 


( 2 ) 




(3) 


V.V+^ = 0, 


(4) 


d(l> 



(5) 
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where 


e = T{p,ipy-e„ 

<P =Oo[Cp(Po/ pV + 8z/do -Cp]. 

Friction is neglected and Q represents dry convective adjustment. The quantities po and 
do are constant and the other symbols have their usual meteorological meaning. 

No along-front variation will be considered here and the two-dimensional domain 
is established in the x-^ plane. The domain is bounded by two horizontal rigid plates and 
periodic boundary conditions are used in x. The top and bottom boundary conditions are 

;=0 at 4 : =0.1, (6) 

and the dependent variables are 


V = u{x, ^, r)i + v(x, ^, Oj, 

(7fl) 

II 

(7b) 

e=e(x,^,t). 

(7c) 


(7d) 


Using the variables in (7) and the condition d/dy = 0, the system (2) - (5) becomes, 
in flux form 
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(8) 


du duu duC d<t> 

dt dx dx ^ ' 


dv duv dvC 


de dud dec ^ ^ 

dt'*' dx d^ 


dud^ 

dx'^d^~^’ 


d; e. 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 


The hydrostatic equation (12) can be vertically integrated and solved for 




(13) 


where 


^0 


and ^ is the value at the lower boundary. The variable 0 cannot be determined from the 
hydrostatic equation because of the rigid upper boundary, since 0 is not known at any 
level. To fix this problem, a dynamic condition is obtained from the continuity 
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equation. Vertically integrating (11) from ^ = 0 to 1 and using (6) yields 

£(«) = 0, (14) 

where the overbar means vertical integration over the whole domain. This implies that the 
mean flow is independent of x. Multiplying (14) by the depth of the layer H establishes 
that the mass flux is constant in space. We further require the mass flux to be independent 
of time to close the system, so that 




(15) 


To use the condition (15) we will rewrite (8) with the use of (13): 


du 

It 


duu du^ 
dx 


^<t>o . 

dx dx ^ 


(16) 


The terms on the right-hand side of (16) can be calculated if V and Q are known, except 
for the unknown term d^dx. Equation (16) is used to predict u without considering 
d^dx and condition (15) is applied to adjust Hu. 

B. MODEL IMPLEMENTATION 

The system of equations is solved with finite differences centered in space and 
time (leapfrog scheme). The nonlinear terms are formulated to conserve the total energy 
(potential plus kinetic). The Arakawa B (Arakawa and Lamb 1977) grid scheme is used in 
the X direction for hydrostatic consistency, so mass and momentum variables are 
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staggered. In the vertical, the energy conserving Lorenz (1960) scheme is applied where 
all variables but vertical velocity are placed in the middle of the layer, to represent the 
average of that layer. The vertical velocity grid points are set between the layers, and 
represent the interface changes. Figure 1 shows the arrangement of the variables in the 
grid system. 



Figure 1. Model grid system following Arakawa and Lamb (1977) and Lorenz (1960). 

The Matsuno scheme is applied periodically at time intervals to prevent solution 
separation. Convective dry adjustment is made when static instability is detected in the 
breakdown of waves. The constants and finite difference increments used in the model 
are displayed in Table 1. The grid spacing in the x direction is variable and is set to give 
the finest resolution in the position of the firontal formation. Figure 2 shows the horizontal 
grid spacing as a function of x. In the region of finest resolution, the horizontal increment 
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is 507 m. The time increment is variable between cases to satisfy the CFL condition 
because the waves propagate at different phase speeds depending on the static stability. 


Table 1. Constants and finite difference increments employed in the model. 


1 Parameter/Constant/Increnient 

Symbol 

Value 

Coriolis parameter 

/ 

10 '^ s'* 

Gravity acceleration 

8 


Reference potential temperature 

Oo 

300 K 

Reference pressure 

Po 

1000 mb 

Vertical domain 

H 

2000 m 

Horizontal domain 

W 


Horizontal space increment (minimum) 

Ax 

507 m 

Vertical space increment 

AC 

1/50 

Time increment 

At 

variable 

Initial static stability 

dGidz 

variable 

Initial temperature anomaly 

AO 

variable 


The numerical solutions are nondimensionalized in the following manner: 


X =xa 


D ’ 


t' = ft , 


6'= 6/AO, 


(u'y)=(u,v)/ 


gAen 


C'=C/ 


gAGH 
\ % J 


where the primes denote the nondimensional variables and is the length scale of the 
initial temperature anomaly AG. The vertical dimension ^ is already in nondimensional 
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Figure 2. H 

resolution (A 
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HI. THEORY 


A mass-momentum imbalance in a statically stable atmosphere will generate 
motion as the fluid tends toward a steady state thermal wind balance (Haltiner and 
Williams 1980; BW96). This geostrophic adjustment happens because the initial 
imbalance excites inertial gravity waves that carry energy away from the region of 
imbalance. The final state is dependent on whether the horizontal length scale ( 0 ^) of the 
region of imbalance is larger or smaller than the Rossby radius of deformation (Lr). The 
final state is determined by the initial mass field for the former case and by the initial 
wind field for the latter. The Rossby radius of deformation is defined as 



where c is the speed of propagation of the appropriate gravity wave (see Haltiner and 
Williams 1980): 


HN 


c = 


K 


and N is the Brunt-Vaisala frequency, defined by: 


which yields 



(17) 
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Unless friction is considered, the gravity waves produce a periodic oscillation about the 
geostrophic balance of constant amplitude in a limited domain. In an infinite domain, the 
gravity waves propagate away from the imbalance region, leaving a state of geostrophic 
balance. In a simple unbalanced initial state, there is a temperature gradient and the 
velocity field is zero. As a result of the imbalance between the mass and momentum 
fields, a transient velocity field develops and helps to further intensify the temperature 
gradient. When the motion is independent of y, u is the divergent part of the wind field 
and V is the rotational part of the wind. The divergent wind causes horizontal convergence 
in the region of maximum temperature gradient which intensifies the gradient and 
enhances upward motion. This nonlinear effect does not happen in quasi-geostrophic and 
linear models. Williams (1972) suggests that the nonlinear results are due to ageostrophic 
effects and these effects tend to accelerate frontal formation. The linear models indicate 
the formation of a temperature discontinuity only as t o® while nonlinearity can 
produce a front within a finite period of time. Another nonlinear effect is the tilting of the 
front toward the cold air, which is important to avoid a superadiabatic lapse rate in the 
warm sector. 

The result of the geostrophic adjustment is different near the upper and lower 
boundaries than away from it. Without the influence of the boundary and for positive 
values of static stability, the geostrophic adjustment is a function of the initial wind field 
(a^ < Lr) or a function of the initial mass field (a^ > Lr). In the former case, the gravity 



waves remove energy from the imbalance region and the temperature gradient decreases 
because the initial wind is zero in the simple case previously mentioned. The larger the 
static stability, the faster the gradient decreases, since the vertical motion plays an 
important role in the temperature change. Near the upper and lower boundaries, the 
vertical motion is much smaller and the temperature changes are much less during 
geostrophic adjustment. When the nonlinear terms are used, the gradient increases near 
the boundaries and the front tilts toward the cold air, in a more realistic fashion. 

Depending on the characteristics of the fluid, the existence of a material boundary 
can cause the temperature gradient to collapse into a front near that boundary, instead of 
oscillate. The formation of a temperature discontinuity in a steady state atmosphere 
depends on the characteristics of the fluid as defined by the Rossby (Ro) and the Froude 
(F) numbers (Williams 1967; BW96). Forf^ 0, they can be defined as 


Ro = 


^gHAe/ep 

/«D 


ylgHAe/d, 

NH 


( 19 ) 

( 20 ) 


Another way of characterizing the flow is by using the Rossby and the Burger number 
(B), which is the ratio Ro/F and can be written 


HN 


/«D ■ 


Equations (18) and (21) can be combined to allow: 


( 21 ) 


15 




( 22 ) 


at n ’ 

and thus B/n Ss, a measure of the adjustment regime for the fluid. 

The definition of a frontal formation in a numerical experiment in not a mere 
objective task. This is particularly true for marginal cases since there is no clear cut 
between the frontal and norifrontal regimes. A useful way to quantify the intensity of the 
front is the d-value {d), used by Williams (1967;1972;1974) and Williams et al. (1992). 
It can be defined as 




AG 


at\de/dx\ 


( 23 ) 


If all the temperature change happened in one grid spacing, as seen in Figure 3, 
d—^ Ax/a'^, as a result of the discretization of the domain. Since J is a measure of the 
temperature gradient and no friction is considered, it is expected to oscillate when no 
frontogenesis occurs and to collapse (d -^Ax/a'^) when a front forms. Since Ax imposes 
a constraint for d in the frontal cases, different grid sizes can be used for the same type of 
flow to assess whether or not a front has occurred. It is also important to verify the 
number of grid spaces involved in a frontal region since the collapse is expected to 
behave as in Figure 3. 

BW96 derived a time-dependent analytic solution for the mass and momentum 
fields that applies to level boundaries where w = 0. They used momentum coordinates 
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-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 


Figure 3. Hypothetical temperature distribution in a frontal zone. The circles represent grid points. 


following Hoskins and Bretherton (1972): 


X = jc + v, 


(24a) 


Z = z, 


(24*) 


(24c) 


where X is conserved (dX/dt = 0), as obtained by substituting (24a) in their y-momentum 
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equation. The use of a domain independent of y allows them to write the momentum 
equations, in nondimensional form, as 


du du du , . 

—+M—+w—-v' = 0, 
at ax dz 

dv' dv dv - 

—+M—+M'^+M = 0, 
at ax ai 


(25) 

(26) 


where v = Vg + v’, and v’ is the ageostrophic velocity. The horizontal pressure gradient 
force is balanced by the geostrophic wind. 

The transformation of (25) into momentum coordinates yields 


du du du , du 

-:rz + u— —hw—-V 

dr dK dz dx 


(dv' dv ^ ^ n 

—+U—+W— 1=0. 
at ox az 


(27) 


Substituting (26) into (27) leads to 


du du , 

—+H’—-V =0. 

dr dz 


(28) 


Following the same process, (26) can be expressed as 


dv' dv ^ 

T—+W—+M = 0. 

dr dz 


(29) 


The solution of v’ for the system (28)-(29), when w = 0 (at 0,1), has the form 


v'=A(X)sin7’ + 5(Z)cosr. 


The initial conditions are « = v = 0, such that on the vertical boundaries 
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v'(X) = -Vg(Z), 


dv'(X) 

dr 


= 0 . 


Therefore, the solution becomes 

v' = -VgCOsr, 

and the total rotational wind solution is 


v = Vg(l-cos7^. (30) 

Blumen and Wu (1994) give the condition for frontogenesis following Hoskins 
and Bretherton (1972) as 


"f-"' 


(31) 


and substituting (30) into (31) leads to the general condition for frontal formation 


(1-C0S7-)|-=1. 


(32) 


For the zero potential vorticity case, the solution for the geostrophic flow can be 
expressed as 


^ j 


(33) 
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Substitution of (33) into (32) gives the particular condition for frontal formation on the 
boundaries when the potential vorticity is zero 

(l-cosr)^ = -2, (34) 

OXq 

since dIdX = 9/8jco- Using an initial temperature distribution of the form 

%M = ]^eTf{RoxQ\ (35) 

equation (34) can be rewritten by going back to dimensional quantities 

(2we)'^ 

= (36) 

Blumen and Wu (1994) used the steady state version of this equation which has no cosine 
to compute the critical value of Ro = 2.03. BW96 included the time dependent part of 
(36) which shows that the minimum Ro is given by 

(/?£>) ^ — —1.435, 

and it occurs at This is the limiting case for large Froude number. 

The authors also derived a more general time-dependent solution for d. The 
dimensional temperature gradient necessary for solving (23) is given, in momentum 
coordinates, by 
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^dd 'l _ A6 ddQ 1 
^dx a^Ro Bxq l-dvjdxQ 


Using (30), (33), and (35) leads to 


(—) = — 
1 _. 




n 


V2 


gABE 

00 


1 'Ncin^ 

--Z 




( 37 ) 




The maximum temperature gradient can be found by differentiating (37) with respect to 
Ro xq and setting the result to zero, which gives 


(Rn^ 1 -2 

w ^7r~nji 




00 


i-z 

J 


sin {fi/2) ^-(Rox,)l 


max 


f 


(38) 


Equation (38) can be solved iteratively for Roxq, which finds the value at the point of 
maximum temperature gradient for a given time. The rf-value is found by substituting 
(38) into (37), and using (23) 


l-2(Roxo)l^ 


(39) 


Equations (38) and (39) also allow the calculation of the time for a front to form. 
Analytical frontogenesis occurs when -> 0 which leads to (i?oxo)^ax =1/2. The time 
for frontal formation can thus be calculated with the use of (38). 
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The time evolution of their solution is shown in Figure 4 for the limiting case. It 
does not consider the barotropic pressure contribution, which is not very significant for 
large temperature forcing associated with frontogenesis. However, this term is important 
for weak forcing. The initial c/-value is 



Figure 4. Time evolution of the 4-value for the analytical solution of BW96. It applies to the limiting case 
of frontal formation, for zero potential vorticity. The critical time in this limit is r = when the analytical 
solution collapses. 

Figure 5 shows some marginal cases of frontal formation as suggested by Blumen 
(personal notes, here denoted Blumen) in a x F plot (Figure 5a is a log-log plot and 
Figure 5b is a linear plot corresponding to the inset in 5a). Blumen extended the analytic 
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theory in BW96 to the stratified solution in Blumen and Wu (1994). The solutions in 
Blumen and Wu (1994) are steady state solutions in geostrophic momentum coordinates. 
The initial state in physical space is the linear geostrophic adjustment solution which 
would result from the initial state used in this thesis. The correction for the time 
dependent factor gives Roj,Fj-Jl. These cases are displayed as circles and the solid 
line is an exponential fit for those points. Frontogenesis occurs for atmospheric flows 
whose Ro and F are located on the concave region of the graph (Region II). The points on 



Figure 5. Blumen’s cases for marginal frontogenesis marked as circles. The solid line is an exponential fit 
for those points. Region I represents the nonfrontal cases and Region II the frontal ones, (a) is a log-log plot 
and (b) is a linear plot corresponding to the inset in (a). 
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the other side of the curve (Region I) represent flows where geostrophic adjustment and 
inertial oscillation occur without collapsing into a front. 
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IV. INITIAL AND BOUNDARY CONDITIONS 


The initial temperature field has an anomaly component which is a function of the 
initial position of the maximum temperature gradient (xq) superimposed on a reference 
temperature that varies with height at a constant positive rate d&dz: 

d= - Z + doix). 

8 

The anomaly component has an error function distribution defined as 

©0 = ^ eTf[ix - Xo )aD \ (40) 

where erf is the error function. The initial position xq is selected for each case in such a 
way as to obtain the minimum i/-value in the region of minimum horizontal grid spacing. 
Lower temperatures lie on the western portion of the domain. Since the amplitude of the 
error function is 2, the amplitude of (40) is given by Aft The temperature distribution 
given by (40) is seen in Figure 6, as a function of x - xq. The initial momentum field is 
zero and thus is out of geostrophic balance with the mass field. 

The assumption that the atmosphere is incompressible and the use of an 
antisymmetric temperature field will lead to antisymmetric solutions on the upper and 
lower boundaries. There is in fact a small difference between the solutions on both 
boundaries due to the different grid spacing such that the fronts will form in regions of 
different resolution. Since the lower front will form in a region of finer resolution, the 
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numerical solution there will be more accurate. Therefore, the lower boundary processes 


will be described in this work, even though it applies to the upper boundary as well. 



Figure 6. Initial temperature anomaly as defined by equation (40), as a function of - j:o- 

The model uses variable grid size in the horizontal to give a fine resolution in the 
region of interest. The horizontal domain is large enough to preserve the region of 
imbalance from lateral boundary effects. The boundaries are far enough so gravity waves 
cannot reach them. Zero vertical velocity condition is used on the vertical boundaries. 
The model uses periodic lateral boundary conditions and two symmetric temperature 
gradient regions to match the temperature in those boundaries. One of the regions is in 





thermal wind balance with the initial wind field and does not change during the 
integration period. Only the unbalanced region is depicted in this study. 
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V. NUMERICAL SOLUTIONS 


The starting point for this investigation is the Blumen’s curve shown in Figure 5. 
The work consists of integrating the model for several cases, whose characteristics are 
described by Ro and F. For convenience, all the parameters of the fluid (listed in Table 1) 
are held constant, except the initial temperature anomaly (A^ and the static stability 
(BOIdz). The nondimensional numbers Ro and F of every case allow the calculation of A6 
and dd/dz, by using equations (19) and (20). For each situation, two different horizontal 
grid spacings are used to help determine whether or not a front occurs. The minimum grid 
spacing in one is twice as large as the minimum for the other. Therefore, the rf-value in 
the former is expected to be roughly twice the value of the latter for a front. If there is no 
front, the values will be the same. 

Figure 7 shows the selected cases (circles) for this study and the curve fit (solid 
line) for Blumen’s points, as displayed in Figure 5. The dashed line represents B = n that 
separates the two different regimes of geostrophic adjustment in large (a^ > Lr) and 
small (a^ < Lr) scale processes, obtained from setting Lr/o'^ = 1 in equation (22). The 
two lines divide the domain in four regions. In regions I and n there is no frontal 
formation, while near-boundary fronts occur in regions El and IV. Regions I and IE are 
characterized by large scale geostrophic adjustment whereas the small scale adjustment 
occurs in regions E and IV. 

Section A presents the analysis of the zero potential vorticity (large F) atmosphere 
to compare with the numerical and analytical solutions of BW96. Examples of the four 


29 



36 CIDO 30 


48 0 000 29 


Figure 7. Distnbution of cases in the present study (circles) and the exponential fit for Blumen’s cases 
(solid line). The dashed line represents B = tc (separation of large- and small-scale geostrophic adjustment). 
The axes are in logarithm scale. The two lines separate the domain in four regions, shown in roman. The 
numbers identify the closest circle and will be refered in the next sections. 

regions shown in Figure 7 are presented in the next four sections. The transition zone 


between large- and small-scale geostrophic adjustment and between frontal and 






nonfrontal regions are discussed in section F. The zero planetary vorticity case is analyzed 
in section G to compare to a large-i?o nonzero case. 

A. ZERO POTENTIAL VORTICITY 

The main characteristic of the zero potential vorticity (F -» «») atmosphere is that 
gravity waves cannot propagate away from the imbalance region due to the absence of 
vertical stability. This means that the mid-level region (away from the upper and lower 
boundaries) does not show a significant change except for the tilting of the isothermals. 
The model levels near the vertical boundaries, however, do suffer significant 
modification. The larger the Ro, the more these levels change. 

Figure 8 shows the time evolution of the rf-value for case 36 {Ro = 1.43). The 
solid line is the numerical value, and the dashed line is the analytic solution obtained 
from equations (38) and (39). There is no front and the temperature gradient oscillates 
with period 27tf(the inertial period). Although the period of oscillation matches, the 
amplitude differs due to the absence of the barotropic term in the analytic solution of the 
d-value. The temperature anomaly and the rotational wind component are displayed in 
Figure 9, for t The divergent component of the wind is zero at that time. 

If the temperature forcing is sufficient to form a front, as in case 38 {Ro = 1.63), 
the inertial oscillation is interrupted because of the temperature jump. Figure 10 depicts 
the time evolution of the d-value for this case. The solid line is the numerical value, and 
the dashed line is the analytic solution obtained from equations (38) and (39). Again, the 
absence of the barotropic term in the analytic solution of the ^f-value causes a difference 
between the numerical and the analytic solutions. For frontal cases, the analytic solution 
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Figure 8, Time evolution of d for the zero potential vorticity case (no front). The solid line is the numerical 
value, and the dashed line is the analytic solution. The period of oscillation is (the inertial period). 




Figure 9, The (a) temperature anomaly and (b) rotational component (v) of the wind for t = (no front). 
The divergent component of the wind is zero at that time. The contour interval in (b) is 0.1 and the zero line 
is removed for clarity. The maximum southerly wind occurs in the lower boundary. 
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precedes the numerical as shown by BW96. The temperature anomaly and the rotational 
wind component, for case 38, are displayed in Figure 11, for t = 2.7/"’. The divergent 
component of the wind tends to zero at that time. 



Figure 10. Time evolution of d for the zero potential vorticity case, with frontal formation. The solid line is 
the numerical value, and the dashed line is the analytic solution for the zero potential vorticity frontal 
formation. 


B. NO FRONT - LARGE-SCALE ADJUSTMENT (REGION I) 

The time evolution of the total temperature gradient of point 48 is presented in 
Figure 12. The fluid initially has Ro = 1.3 and F - 15.3. The temperature gradient 
increases at the boundaries as the zone of maximum gradient moves toward the warmer 


33 






Figure 11. The (a) temperature anomaly and (b) rotational component (v) of the wind for t = 2.7/'*, with 
frontal formation. The divergent component of the wind tends to zero at that time. The contour interval in 
(b) is 0.1 and the zero line is removed for clarity. The maximum southerly wind occurs in the lower 
boundary. 

air in the lower boundary and toward the colder air in the upper one (the axis of 
maximum gradient tilts so that the colder air is underneath). The oscillatory nature of the 
temperature distribution can be noticed by comparing the cross-sections at f = 0 and 
t = This case is very similar to that shown in section A, despite the large difference 
inF. 

Figure 13 displays the value evolution with time for the levels (a) ^=0 and (b) 
^=0.5. The changes in the maximum temperature gradient at the former level are much 
larger than those in the latter. The static stability is very small and gravity waves are very 
slow, which results in very little energy propagating away from the initial imbalance 
region. The main cause of temperature distribution change is due to the divergent part of 
the wind (u), near the vertical boundaries. Figures 14 and 15 show the wind distribution 
(m and V respectively) associated with the temperature distribution in Figure 12. The 
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Figure 12. Time evolution of the total temperature cross-section for point 48 in Figure 7, for (a) r = 0, (b) 
t = l.Sf\ (c) t = nf\ (d) t = 4.7/', (e) t = 2n/-', and (f) t = 9.8/-'. 
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Figure 13. Time distribution of value (case 48) for levels (a) ^ = 0 and (b) 0.5. 


M-component works both to intensify the temperature gradient when it is small and to 
diminish it when it is strong, causing the oscillation seen in Figure 13. This component 
shifts sign in the points of maximum and minimum </-value. The rotational component 
anomaly (v), on the other hand, attains its maximum (minimum) when and where the 
temperature gradient is maximum (minimum). Both components are nearly zero in the 
vicinity of the level ^=0.5. This relationship is depicted in Figure 16 that shows the time 
series for the point where the </-value is calculated (maximum temperature gradient), in 
the lower boundary. As a result of the inhomogeneous surface convergence, the 
temperature distribution near the boundaries loses the antisymmetric characteristic as 
seen in Figure 17, at r = 0 and t = nf^. 

The use of a grid spacing with the minimum Ax twice as large as the previous one 
results in similar fields and cf-value. The other points in this region also show analogous 
behavior except that the points with larger Ro have a small decrease in the temperature 
gradient near 0.5. In these cases, the oscillation can be noticed in the tilt of the 
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Figure 14. Time evolution of u for point 48, 
for (a) r= 1.8^^*, (b) f = 4.7/-and ic)t = 9.Sf\ 
The M field is nearly zero for the other times in 
Figure 11, The isoline interval is 0.1. In the 
lower boundary, u is positive in (a), and negative 
in (b) and (c). The zero isoline is removed for 


Figure 15. Time evolution of v for point 48, 
for (a) t = 1.8/'\ (b) t = and (c) t = 4.7/‘\ 
The V field is zero for f = 0 and t = and 
very similar to (b) for t = 9.8/ In the lower 
boundary, v is negative in the three frames. The 
zero isoline is removed for clarity. 













isotherms. Another difference is the change in the period of oscillation that decreases as 
Ro increases, but no difference is observed when F changes. 





Figure 16. Time series for (a) </-value, (b) u, and (c) v in the point of maximum temperature gradient, in the 
lower boundary (case 48). 


C. NO FRONT - SMALL-SCALE ADJUSTMENT (REGION II) 

Figures 18 and 19 show the time variation of the temperature field of point 33 (see 
Figure 7), which has Ro = 2.975 and F = 0.602. The former presents the total temperature 
whereas the latter depicts the temperature anomaly field. The main feature of this case is 
the tendency of the temperature gradient, away from the boundaries, to be strongly 
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Figure 17. Temperature profile for the lowest level, at f = 0 (solid line) and t = 7^* (dashed line). 

damped. This is expected since in this region the steady state solution is a function of the 
initial momentum field, which is initially zero. Therefore, in the long run, the temperature 
gradient will tend to vanish. 

The temperature gradient in the lower boundary, on the other hand, behaves the 
same way as in region I, where the key player is the divergent wind field. Figure 20 shows 
the d-value for the lower boundary, where the characteristic oscillation is seen. The 
period is much smaller than that of region I. The value for the mid-level region is not a 
good indication in this case because the gravity waves become very fast and distort the 
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Figure 20. Time distribution of d-value at the lower boundary for case 33. 


isotherms considerably. This distortion can mask the larger scale temperature gradient in 
the value calculation. 

Another important difference is that the momentum field lags the mass field. 
Although the spatial wind distribution in this case (not shown) is similar to those in 
region I, the lagging is substantial. This is due to the rapid change in the temperature 
field by the gravity waves compared to the inertial period (27^'^). Figure 21 shows the 
time series for (a) d-value, (b) u, and (c) v in the point of maximum temperature gradient, 
in the lower boundary. 
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Figure 21. Time series for (a) rf-value, (b) u, and (c) v in the point of maximum temperature gradient, in the 
lower boundary (point 33). 

For points with larger Ro the period of oscillation decreases considerably 
compared to The propagation speed of the fastest gravity waves increases and 

therefore the geostrophic adjustment process becomes much faster. As a result, the 
temperature gradient away from the vertical boundaries decreases faster. 

D. FRONTOGENESIS WITH LARGE-SCALE ADJUSTMENT (REGION IE) 

The distinct characteristic for this region is the temperature collapse into a frontal 
zone in the upper and lower boundaries, while little change is observed away from the 
boundaries. The temperature anomaly at t = 2.7/'\ in the lowest level, for case 29 
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Figure 22. Temperature anomaly in the lower boundary for case 29 at r = 2.7/'*. The frontal zone is the 
steep slope located at jc - Jto = 1.5 Cp. 

(Ro = 1.630 and F = 15.302), is depicted in Figure 22. The large temperature change in a 
very narrow zone characterizes the collapse, which causes the model solution to be 
discontinued. It can be noticed that only a fraction of the total temperature anomaly is 
present in the frontal region (at a: - = 1-5 a'^), while the remainder goes into a very 

smooth gradient. The frontal zone near the surface moves rapidly toward the warmer air, 
causing the front to tilt toward the colder air. Figure 23 shows the J-value changes at the 
surface for case 29, which agrees very well with BW96. The dashed line has the 
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minimum grid size twice as large as the solid one. The grid size imposes a constraint on 
the d-value when the collapse occurs. Therefore, the minimum rf-value in the first case is 
roughly twice that of the latter. The ^/-value at level ^ = 0.5 (not shown) indicates very 
small changes in the temperature gradient at this level. 



Figure 23. Time variation of the d-value for the lower boundary, for case 29. The minimum d-value is 
lower for the smaller grid scheme (solid line) than for the larger (dashed line). Note the fast decrease and 
the collapse near t = 2.7/■^ 

Figure 24 shows the evolution of the temperature anomaly field for (a) f = 0, 
(b) r = 1.1 /■*, (c) r = 1.8/'*, and (d) t = 2.7/'*. The momentum field for the same time is 
shown in Figures 25 (m component) and 26 (v component). This result is similar to the 
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Figure 24. Time evolution of the temperature anomaly cross-section for point 29 in Figure 7, for (a) t = 0, 
(b) r = 1.1 f *, (c) t = 1.8/-*, and (d) t = 2.1 f\ 


case showed in BW96, with The same relationship between mass and momentum 

fields observed in region I also occurs in this region. The rotational component increases 
in phase with the gradient increase, due to the Coriolis force. The divergent component, 
on the other hand, increases initially because of the unbalanced pressure gradient, but is 
very small when the temperature collapses. This is more evident in Figure 27, which 
displays the time series of d-value, u and v. 
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Figure 25. Time evolution of u for point 29 in Figure 26. Time evolution of v for point 29 in 

Figure?, for (a) t = 1.1 /■’, (b)r=1.8/-\ Figure 7, for (a) t = 1.1 /'*, (b) r= 1.8/-*, 

and (c) t = 2.7 /■*. The zero isoline is removed and (c) t = 2.7 /'*. The zero isoline is removed 


for clarity. 


for clarity. 














Figure 27. Time series for (a) rf-value, (b) u, and (c) v in the point of maximum temperature gradient, in the 
lower boundary for case 29. 

E. FRONTOGENESIS WITH SMALL-SCALE ADJUSTMENT (REGION IV) 
As in region DI, the temperature gradient collapses into a front but more rapidly. 
Figure 28 presents the t/-value time variation for case 44 (Ro = 4.000 and F = 0.701), 
which is not different from the cases in region HI, except for the time scale. Again, the 
dashed line represents the case where the grid spacing is twice the one in the solid line. 
Figure 29 shows the temperature anomaly field in the lowest level, close to the lower 
boundary, at r = 1.25/'^ There is only part of the temperature change involved in the 
frontal zone as before. The main difference between this region and the previous one is 


48 





the strong damping in the mid-level temperature gradient as depicted in Figures 30 
and 31. 



Figure 28. Time variation of the <f-value for the lower boundary, for case 44. The minimum d-value is 
lower for the smaller grid scheme (solid line) than for the larger (dashed line). The frontal formation occurs 
near r= 1.25/'*. 


The u and v fields are shown in Figures 32 and 33. Note that unlike the case in 
region IE, the divergent wind is not negligible at the time of the frontal formation. There 
is a lag in the momentum field with respect to the mass field, as seen in Figure 34. This 
time series show the correlation and the lag between those fields. 
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Figure 29. Temperature anomaly in the lower boundary for case 44 at t = 1.25 The frontal zone is 
located j: - = 2.2 Cp. 

The other cases in this region present the same behavior as case 44. A difference 
is the time scale involved in the frontal formation in the lower boundary. The other 
distinct characteristic is the geostrophic adjustment process in the midlevel zone. This 
adjustment is stronger for larger Ro. 

F. TRANSITION ZONES 

The transition zone between the large- and small-scale geostrophic adjustment 
regions is not very large. Although the cases in each region show some variation in terms 
of geostrophic adjustment as they approach the transition zone, this change is much more 
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Figure 30. Time evolution of the total temperature field for point 44 in Figure 7, for (a) f = 0, 
(b) t = 0.63 f\ (c) t = 0.90 f\ and (d) t = 1.35/-*. 

abrupt near the line B = n. This is true both when changing between regions I and n, and 
between DI and IV. Cases 2 and 6 (see Figure 7) are located in both sides of the line 
B = n. They show very distinct characteristics in the mid-level temperature gradient. The 
former lies slightly to the right of the curve (small-scale region), and undergoes a large 
change in the temperature gradient as gravity waves move away from the imbalance zone. 
This changes are very similar to those depicted in Figures 18 and 19. The latter case is 
located in the region of large-scale geostrophic adjustment, and little change occurs in the 
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temperature gradient away from the boundaries. This behavior is analogous to the case 
showed in Figure 12. 


Figure 31. Same as in Figure 30, except for the temperature anomaly. 


The transition from nonfrontal to frontal regions is not clear cut. The approach 
here is to consider that frontogenesis has occurred if all the situations below have 


occurred in the lowest level of the model (near the lower boundary): 
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Figure 32. Tune evolution of u for point 44, for 
(a) t = 0.63f\ (b) t = 0.90/-', and (c) t = 1.35/-‘. 
The zero isoline is removed for clarity. 


Figure 33. Time evolution of v for point 44, for 
(a) t = 0.63/-', (b) t = 0.90f'\ and (c) t = 1.35/-‘. 
The zero isoline is removed for clarity. 













•^/-value at critical time is bounded by the grid size; 

•rf-value tends to zero at the critical time; 

• there is no oscillation in the temperature gradient after the critical time; and 

• there is a large temperature change in a few grid points. 


d 





tf 


Figure 34. Time series for (a) rf-value, (b) u, and (c) v in the point of maximum temperature gradient, in the 
lower boundary for case 44. 


The first criteria is evaluated by comparing the J-value for a given characteristic 
flow (Ro, F), run with two different grid spacings. Figure 35 presents the </-value for 
Ax = 0.01 (dashed) and Ax = 0.005 (solid), for case 20. The ratio of the two Ax is 1:2 
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while the ratio of the rf-values is 1:1.7. This is very close considering that finite difference 
models have different precision for different Ax and that this is a non-linear model. 
However, not all the change in temperature is in the frontal zone and thus the ratio 
between the different Ax is not a precise measurement of the ratio between the two 
J-values. Nine cases were found to be in this transition zone. The average ratio of J-value 
is 1:1.6 with a standard deviation of 6%. 



Figure 35. Time evolution of J-value for case 20 for Ax = 0.01 (dashed) and Ax = 0.005 (solid). The front 
collapses near t = Tcf’. 
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Model runs have shown that d = 0.10 marks the beginning of the frontal region. 
There is no frontal formation when > 0.10 and a front has formed when d < 0.10. 
Analysis of case 20 {Ro = 1.543 and F = 7.999) shows a good consistency with the above 
characteristics. It is located between a nonfrontal and a frontal case and has a minimum 
d-0.\0 which is much smaller than the initial J-value. It is important to notice that this 
limiting d = 20Ax and for stronger frontal cases d—^ Ax. 

The frontal formation causes the collapse of the temperature gradient and thus 
interrupts the oscillation observed in nonfrontal cases, which leads to the third criteria. It 
can be seen in Figure 35 that the temperature gradient in case 20 does not oscillate as in a 
nonfrontal situation. After the temperature collapse, there is no meaningful numerical 
solution and the model output is discontinued. The last criteria consists in observing how 
close the temperature distribution resembles the ideal case shown in Figure 3. The 
temperature distribution in the lowest level of the model is shown in Figure 36 (case 20). 
The circles represent the grid points and indicate the number of grid spaces that are 
involved in the temperature jump. This jump usually happens in three to four grid spaces, 
which is a rather small frontal zone with horizontal scale of 0(0.020^) for the smaller 
grid spacing and 0(0.04^0) for the larger one. It is important to notice that the scale 
continues to fall as long as Ax is decreased. 

An empirical exponential-type curve that establishes the limit for frontal 
formation is presented in Figure 37. This limit is defined by the following relation 
between Ro and F: 

F = ae~^^^^^^\e, (41) 



where a = 25.0, b = 3.90, c = 1.53, d = 0.23, and e = 0.37. Frontogenesis occurs for F 
larger than the RHS of (41). An important characteristic is the asymptotic convergence to 
Ro = 1.53. The lower part of the curve shows the tendency for asymptotic convergence to 
a limiting F. To evaluate the convergence for larger Ro it would be necessary to use a 
finer resolution model with larger domain. The stability increases excessively and causes 
the gravity waves to propagate very fast. As a result, the front also moves very fast and 
requires a larger zone with finer resolution. The domain must be increased so that the 
gravity waves will not reach the lateral boundaries. The curve seems to converge to a 
value close to F = 0.37. 



Figure 36, Temperature distribution in the lowest level of the model for case 20. The circles represent the 
grid points. The jump occurs in few grid points. 
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G. ZERO PLANETARY VORTICITY 


The planetary vorticity (f) is very important for small Ro as seen in sections B 
and D. The period of oscillation of the temperature gradient is nearly inertial (27^'^). As 
Ro increases, however, the oscillation departs from the inertial value, as in sections C 
and E. For large values of Ro, the planetary vorticity plays a very small role. Case 57 (see 
Figure 7) has = 9 and F = 0.5, and was integrated using/= lO'^s'*. It is compared to a 
similar case with/= 10"^ s’* (Ro oo). Figure 38 displays the time evolution of the d- 
value for/= 10’^ s’* (solid line), and /= 10’* (dashed line), where the horizontal axis is 



Figure 38. £?-value evolution for /= 10"^ s ' (solid line), and /= 10 * s'' (dashed line). The former slightly 
lags the latter. 
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the dimensional time (in hours). The two curves start at the same point and slightly 
separate, beginning at r = 0.3 hs. Although the absence of planetary vorticity causes the 
front to form a little earlier than the nonzero case, the 6 and u fields behave similarly in 
both situations. The two cases are different because Ro is not large enough in the former 
case and Ro —>oo as/—> 0. This case is also treated for F -»<» in studies by Simpson and 
Linden (1989) and Kay (1992). 

The temperature and divergent wind component fields for this two cases, close to 
the time of frontal formation (not displayed), show very subtle differences. The only 
significant difference for large Ro atmosphere is in the rotational wind component. There 
is no thermal wind balance due to the lack of planetary vorticity. Therefore, the rotational 
part of the wind is zero when/—> 0. 
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VI. DISCUSSION 


The geostrophic adjustment using the present model is very consistent, and it 
happens in the levels where the influence of the upper and lower boundaries is small. The 
initial unbalanced state excites gravity waves that propagate away from the imbalance 
region, for a very large domain. The mass and momentum fields tend to a geostrophic 
balance as the divergent field vanishes. During the geostrophic adjustment process, the 
temperature and the wind fields oscillate about the steady state solution, as the gravity 
waves propagate away from the imbalance region. 

A limit between the large-scale > ^r) small-scale < Lr) geostrophic 
adjustment regimes is set by fi = it. Therefore B>Tt leads to an adjustment that depends 
on the initial wind field (small-scale process), and B <n depends on the initial mass field 
(large-scale process). The transition zone between the two regimes is small, as discussed 
in chapter V, section F. In the type of initial state here presented, large-scale regime tends 
to have small changes in the temperature gradient in the levels distant from the vertical 
boundaries, except for the tilting of the isothermals. When F -» oo, there is just the 
inertial oscillation. In the small-scale regime, on the other hand, the fields tend to a state 
of no temperature gradient and no motion, because the initial state is at rest. If an initial 
imbalance consisted of a uniform horizontal temperature field and a nonzero wind field, 
the geostrophic adjustment would lead to a different final state. The large-scale final state 
would be zero and the small-scale final state would consist of a temperature field in 
thermal wind balance with the rotational wind field. 
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The temperature gradients in the lower and upper boundaries are not driven by 
the Rossby radius of deformation (Lr). The temperature behavior near the vertical 
boundaries is independent of Bln. Rather, the temperature change depends on the 
intensity of the divergent part of the wind. Convergence near the maximum temperature 
gradient is required for further intensification. The intensity of the divergent wind is a 
function of the initial state of the atmosphere, as defined by the pair {Ro, F). A proposed 
limit for frontal formation is defined by the curve presented in Figure 37. It differs from 
Blumen’s curve mainly for large Ro. Figure 39 displays both curves in a log-log plot. The 
curve suggested in this work (solid line) seems to converge to a limiting F = 0.37. 
Blumen’s curve (dashed line) suggests the convergence would occur to a smaller Froude 
number and not so rapidly. 

Near the upper and lower boundaries, the temperature gradient also oscillates 
when a front does not occur. The period of oscillation (7) is nearly inertial. For points 
distant from the limiting curve, with no frontal formation, the model is nearly linear and 
the £?-value oscillates smoothly with the fundamental period of oscillation. For points 
closer to the curve, however, the nonlinear effects cause high frequency oscillations 
superimposed on the fundamental frequency and the d-value is no longer smooth. 
Figure 40 presents the time evolution of the temperature gradient for cases 1 (Ro = 2.215, 
F = 0.685) and 2 (Ro = 2.476, F = 0.765), the latter being closer to the limiting curve. 
Case 1 oscillates smoothly while case 2 contains high frequency oscillations. Figure 41 
shows the calculated rf-value variation for case 36 (large F), as well as the distribution 
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Figure 40. Time evolution of the temperature gradient value) for cases (a) 1, and (b) 2. Case 2 is 
significantly modified by the nonlinear terms. 

The fundamental period of oscillation is not constant for all situations, and seems 
to depend on Ro, F and N. For Ro < 1.53 (the limiting Rossby number), the oscillation is 
inertial and The period of oscillation changes for Ro > 1.53, as shown in 

Figure 42, where the horizontal axis is RolF and the vertical is NT. The curve seems to fit 
all the cases except the ones located below the limiting F (not shown), because of the 
very high stability. The half-period of oscillation was assumed to be the time the d-value 
takes to achieve its first minimum, starting at the origin. 
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Figure 41. d-value variation for case 36 with (a) regular calculation, and (b) after filtering frequencies 
higher than the inertial frequency, using a Chebychev filter. 

When the pair (Ro, F) is located over or in the concave side of the limiting curve, 
the temperature gradient no longer oscillates. A jump discontinuity forms which 
characterizes frontogenesis. The time ( 4 ) for a front to form in a zero potential vorticity 
atmosphere can be calculated with the help of equations (35) and (36). These equations 
do not consider the barotropic pressure gradient and were derived for the zero potential 
vorticity case. The barotropic term is not significant for strong forcing associated with 
frontal formation and becomes less important form stronger cases (BW96). As a result, 
larger Ro fluids tend to have smaller difference between the analytic (ta) and the 
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Figure 42. Variation of the oscillation period with Ro, F, and N, for Ro 1. 1.53, and F > 0.37. These cases 
(circles) do not produce fronts and oscillates with period T. The solid line is a cubic spline interpolation for 
those points. 

actual time (tc) for a front to form. However, for nonzero potential vorticity situations, 
this relationship is not so straightforward. Stronger fronts form much more rapidly than 
weaker ones. In this study, the time frontogenesis occurs is considered the intercept of the 
extrapolation of the almost linear part of the if-value curve with the time axis, as 
presented in Figure 43 (for case 22). This procedure is more consistent with the analytical 
solution, as presented in Figure 4, and it avoids the change in the slope of d caused by 
using a discretized domain. The ratio tjtc decreases when F decreases but also decreases 
for weaker frontal cases. This means that the ratio is larger for points located far from the 
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Figure 43. The time for a front to occur is obtained by extrapolating the almost linear part of the <f-value 
curve until it intercepts the time axis. In this case, a front is considered to form at t = 2.39/'*. 


limiting curve (in the concave region). Figure 44 presents the ratio tjtc as a function of 
F* for the cases analyzed in this work. The crosses represent the cases with S<d< 12 
(weak ffontogenesis) and the circles, d < S (strong ffontogenesis). The solid line is a 
second order polynomial fit for the weak cases. The points along the limiting (Ro, F), 
marked with crosses, put a lower bound for the tJtc ratio. Stronger frontogenetical cases 
depart from this lower limit as the barotropic pressure gradient becomes less important. 

Another important characteristic of the frontal temperature distribution is the 
relatively small amount of the total temperature change involved in the jump. The 
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Figure 44. The ratio tjtc for the limiting frontal cases (crosses) establishes a lower bound for that ratio. 
The solid line is a second order polynomial fit for these cases. Stronger frontogenetical cases (circles) tend 
to lie above the limit. 


temperature gradient near the upper and lower boundaries change due to the divergent 
wind. The point of maximum gradient is initially located in the center of the imbalance 
region and shifts toward the warmer air as shown in Figure 29. The temperature field is 
modified in such a way that part of the temperature gradient becomes steeper and 
collapses into a front. The remaining part of the gradient becomes smoother. The analysis 
of all cases where frontogenesis occurred show no trend in the data. The average of the 
amount of temperature involved in the jump (normalized by is 0.175 with a standard 
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deviation of 0.045 (for a sample of 28 fronts). This amount does not seem to be changed 
for different parameters and characteristics. The variability is probably explained by the 
fact that model outputs are recorded at an interval At = 0.09/"V It means that the front 
may be formed at a different time than the recorded time. A source of error is the 
difficulty of determining where the jump starts and finishes. The criteria used here was 
the use of the largest changes in the slope of the temperature distribution. 

The change in the position of maximum temperature gradient is also not constant 
for all frontal cases. The speed of the point of maximum temperature from the initial 
position until the jump forms is presented in Figure 45, as a function of Ro. The crosses 



Figure 45. Speed of propagation (Vp) of the position of maximum temperature gradient for frontal cases 
(crosses) as a function of Ro. The solid line is the first order polynomial fit for the crosses. 
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represent the frontal cases in this work and the solid line is the first order polynomial fit 
for those points. The larger the Rossby number, the faster the movement of the 
temperature gradient. The speed is calculated using the change in the position of the 
if-value and the time for frontal formation (tc). 
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VII. CONCLUSIONS 


The numerical model of Williams et al. (1992) is used to study frontogenesis from 
unbalanced initial conditions. The dependent variables are assumed to be independent of 
y. The hydrostatic Boussinesq primitive equations are used with no diffusion of heat or 
momentum. The model is bounded at the top and bottom by rigid planes. Periodic 
boundary conditions are used in the horizontal. Stretching horizontal coordinates are used 
to give fine resolution in the region of frontal formation. The lateral boundaries are placed 
far enough from the imbalance region to avoid wave reflection. The initial temperature 
field has a constant vertical stratification. 

The initial imbalance is obtained by allowing a horizontal temperature gradient to 
exist while the initial wind is zero. In a stably stratified atmosphere, gravity waves are 
excited and propagate away from the imbalance region, provided no reflection occurs at 
the lateral boundaries. Therefore, the atmosphere tends toward a geostrophic balance 
away from vertical boundaries. Near these boundaries, the temperature gradient oscillates 
or it collapses into a front, depending on the initial Rossby (Ro) and Froude (F) numbers. 

A relationship between Ro and F is established to limit situations where a front 
may or may not form. Therefore, if the initial Ro and F of an imbalance region can be 
determined, the unbalanced contribution for future development of discontinuities can be 
estimated. This relationship also shows how intense the process is, depending on the 
distance of the point (Ro, F) from the limiting curve. The scale of the geostrophic 
adjustment far from the boundaries is likewise determined by Ro and F. Large scale 
adjustment happens when F > Ro/n and the fluid adjusts to the initial mass field. The 
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adjustment depends on the initial momentum field when F < Ro/it, and the initial 
imbalance scale is small compared to Lr. 

Large Ro and small F fluids are more divergent and less dependent on the inertial 
period, which leads to faster firontogenesis and more rapid geostrophic adjustment. A zero 
planetary vorticity (Ro —> <») case is compared to a large Ro flow with/= 10*^ s'* (the 
other parameters held constant), and the results show very subtle differences. The period 
of oscillation for the nonfrontal cases is inertial when Ro < 1.53, and depends on Ro, F, 
and N for larger Ro. The amount of temperature involved in the frontal zone is about 17 
% of the temperature anomaly. This differs from the simple inviscid Burger’s equation, 
where almost all the temperature collapses into a front. 

For large F Blumen’s curve is very close to the critical curve in our study, since 
his value for large F is = 1.45. This is also shown by BW96 who are considering the 
zero potential vorticity case. In this case the critical value should be veiy close because 
the initial conditions are the same. The latter are the same because the initial state is not 
changed by geostrophic adjustment. The difference between the Blumen’s curve and our 
curve for large F is caused by the barotropic pressure contribution which is not included 
in the analytic theory. As F decreases the Blumen’s curve falls below and to the left of 
our curve. Although the Blumen’s curve is based on a weaker initial temperature field (it 
is the same along the upper and lower boundaries) this does not seem to be important 
because his curve is still to the left of our curve. If this effect were more important it 
would require larger Ro and a larger F. Apparently the geostrophic adjustment in the 
numerical model brings the temperature close to the initial field used by Blumen. 
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The larger difference between the Blumen’s curve and our curve in the lower part of the 
Ro-F diagram must be caused by the barotropic pressure gradient term. 
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